Thermodynamic Theory of Phase Transitions in Driven Lattice Gases 
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We formulate an approximate thermodynamic theory of the phase transition in driven lattice gases 
with attractive nearest-neighbor interactions. We construct the van der Waals equation of state for 
a driven system where a nonequilibrium chemical potential can be expressed as a function of density 
and driving field. A Maxwell's construction leads to the phase transition from a homogeneous fluid 
phase to the coexisting phases of gas and liquid. 
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Introduction. - Macroscopic properties of systems in 
equilibrium are described by thermodynamic potentials, 
like entropy or free energy, which can be derived from 
the microscopic properties through the Boltzmann dis- 
tribution. The ultimate triumph of this formalism lies 
in describing phase transitions, arguably the most inter- 
esting phenomena known to occur in various interacting 
many-particle systems. 

A phase transition can also occur in a system with 
a nonequilibrium steady state (NESS) which exhibits a 
steady current. However, driven systems have so far re- 
sisted attempts to construct a general formalism similar 
to that in equilibrium |l| . Understandably, studies in this 
direction have focused on getting insights from simple 
model-systems [2j-|8[. One such model for systems hav- 
ing a NESS is the driven lattice gas (DLG) |9| which has 
become a paradigm in nonequilibrium statistical physics, 
analogous to the paradigmatic Ising model or equilibrium 
lattice gas (ELG). Although the DLG has been studied 
extensively in the last couple of decades and the various 
properties concerning the nonequilibrium phase transi- 
tions are fairly well known [lfj, lll| . a thermodynamic 



theory is still missing even for this one of the simplest 
models of driven interacting many-particle systems. 

In this paper, we formulate an approximate thermo- 
dynamic theory which not only captures various macro- 
scopic properties but also describes the phase transition 
in the driven lattice gases with attractive interactions. 
We construct a mean-field (MF) van der Waals equation 
of state for a driven system with a chemical potential 
/i(n, E) expressed as a function of density n and driving 
field E. The quantity /i is identified using the concept 
of equalization of an intensive variable, like equilibrium 
chemical potential, for a driven system kept in contact 
with the corresponding non-driven one. Then we use 
the Maxwell's construction, familiar in equilibrium for 
constructing a concave (or convex) thermodynamic po- 
tential to describe the phase-coexistence, to explain the 
phase transition from a homogeneous fluid phase to the 
coexisting phases of gas and liquid. Our theory is in re- 
markable agreement with the numerical observations. 

Model. - We consider a model, introduced earlier in 
[l2j, of two lattice gases, one driven with volume V\ and 



the other non-driven with volume V2, exchanging parti- 
cles through a small contact at V\ and V2, respectively. 
The energy H of the two systems combined is given by 

H = K 1 J2v{ r i) r l( r i') + ^2Ef( I '2)')( r 2') where sums 
arc over nearest-neighbor sites with ri,ri' £ V\ and 
T2,r2' £ V2, K\ and K^ the interaction strengths of the 
pair-potentials among particles in systems 1 and 2, re- 
spectively, and ry(r) the occupation variable taking val- 
ues only 1 or given the site r is occupied or unoccu- 
pied, respectively. We choose the jump rate w(C'\C) 
from a configuration C to C according to the local de- 
tailed balance condition [9(: the jump rate for a parti- 
cle from a site r to its unoccupied nearest neighbor r' 
obeys w(C'\C) = w{C\C) exp[-Aff + E{x' - x)} where 
AH = H(C') - H(C), E is the driving field along the 
x-direction, and x and x' are ^-components of r and 
r' (krjT — 1, ks the Boltzmann constant, T tempera- 
ture). We consider two-dimensional systems with peri- 
odic boundaries in both directions and choose E = E\ 
when r,r' £ V\, and E = otherwise. Also, we con- 
fine ourselves to the cases where the combined system is 
particle-hole symmetric with K\ = K2 = K and consists 
of particles with attractive interaction of strength K < 0. 
For Ei = 0, the combined system, an equilibrium lat- 
tice gas, has the Boltzmann distribution. However, for 
Ei 7^ 0, there is a current in system 1 in the steady state 
with a probability distribution unknown in general. 

Mean-field theory. - Defining the quantity w a / a as the 
conditional average of the jump-rate of a particle from 
system a to system a' if a contact site in system a is 
occupied and the corresponding contact site in system a' 
is unoccupied with a, a' = 1,2 and a/a', we get 
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where n a is density at the contact site in system a. We 
define the conditional average, which will be useful later, 
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where T) a = J2i=i Vai IS tne sum over variable 77™", the 
occupation variable at the zth nearest-neighbor site to 
the contact site in system a, 77 = or 1, and P{fj a \r]a ) 



is the conditional probability of fj a given a fixed value of 

(c) 

i]a , the occupation variable at the contact site in sys- 
tem a. Now we use a mean-field (MF) approximation 
for the conditional jump-rate w a > a ~ exp[—(K a i(fi a i)o — 
K a (fj a )i)/2] where the quantity in the round brackets in 
the exponent is the conditional average of difference in 
energy between the final and the initial configurations, 
given that a particle jumps from system a to system a' . 
Therefore we get from Eq. [TJ 

n 2 c Jif 3 (rtJ,+(j»J n )/2 = "i e Jfi((^> 1 +<^> )/a - 

1 — 712 1 — Til 

(3) 
Here we have implicitly assumed that correlations be- 
tween the two systems across the contact are negligibly 

(c) 

small and consequently n a ss n a where n a is the bulk- 
density in system a. Finally, putting K\ = Ki = K 
and then taking logarithm, Eq. [3] can be rewritten in the 
more illuminating form, 
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which constitutes the basis of the following analysis. Now 
one can readily identify the l.h.s of Eq. [4] as the chemical 
potential /x 2 = ln[n 2 /(l - "2)] + K ((fj 2 ) 1 + (fj 2 ) ) /2 of 
an equilibrium system, in this MF approximation, with 
density n 2 and interaction strength K. 

Recent studies of DLGs have revealed a simple ther- 
modynamic structure where, in a large parameter space 
and to a good approximation, one could define an inten- 
sive variable, like equilibrium chemical potential, which 
equalizes upon contact 12j. At this point, we use this 
concept of assigning the chemical potential /j of the equi- 
librium system 2, on the mean-field level, to the driven 
system 1. We verify it by checking the fluctuation- 
response relation d{Ni)/dfi = <r 2 N where d(N 1 )/d/j, is 
the compressibility and a 2 N = ((iVf ) — (^i) 2 ) is the fluc- 
tuations in particle-number i\q of driven system 1 which 
is in contact with a particle reservoir equilibrium system 
2. This fluctuation relation is indeed remarkably well 
satisfied as seen in the top panel of Fig. [1] 

van der Waals equation of state. - The above obser- 
vation leads us to identify the r.h.s of Eq. 0J on the 
mean-field level, as the chemical potential /z(n, E) for a 
driven system with density n and driving field E, i.e., 
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where we drop the subscript of the occupation variable i). 
The dependence of fj, on E enters through the conditional 
average of f) which can be written explicitly in terms of 
the nearest-neighbor correlation u = (r](r)ij(r')), with r 
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FIG. 1: Top panel: Simulation results for a 20 x 20 system 1 in 
contact (in the 2x2 contact region) with a 250 x 250 equilib- 
rium system 2, a particle reservoir. Integrated compressibility 
I x = f M (<9(iVi)/<9//)d/i' (squares) and integrated fluctuation 
I a = J M ((j% )dfi' (circles) as a function of chemical potential 
fj, for a driven system (red points) with K — — 1, E — 6 as 
well as an equilibrium one (blue points) with K — — 1 where 
the systems are in contact with equilibrium reservoirs with 
the same respective K. Bottom panel: The MF results where 
driven and the corresponding equilibrium case are compared 
(see main text after Eq. I10|l . All curves obtained in the fluid 
phase. 



and r' two nearest-neighbor sites, 
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where P(rj^) is the probability of the occupation vari- 
able ?/ c ) at the contact site. Now using P(jf c ^ = 1) ~ n, 
P(v = 0) w (1 — n), (rf c 'r)) ~ 4w + where u + = (u|| + 
uj_)/2 the average nearest- neighbor correlations in the 
bulk with My and u± the nearest-neighbor correlations re- 
spectively along and perpendicular to the direction of the 
driving field E, and finally defining the average nearest- 
neighbor correlation function c + (n,E) = (u+ — ?i 2 ), we 
obtain from Eq. [51 
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where the function M(n, E) is defined as 
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The quantities uu, u± and u + depend on E. Note that, 
in deriving Eqs. [7] and [51 we have used that the den- 
sity, and the quantities u\\ and u±, arc approximately 
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FIG. 2: Nearest-neighbor correlation function c+(n, E) for 
DLG and ELG compared as a function of density n for 
K = -1 and if = -0.5. Note that c+(n,S) < c+(n,0). 
The magenta and sky-blue lines are plotted, for if = —0.5 
and K = — 1, using the functional form of Eq. [TO] Inset: 
Scaled c+(n, 0) and C-|-(n, .E) are reasonably well collapsed 
on each other by using the approximate form of c+ (n, _B) ~ 
A(E)c+(n, 0) where j4(_B = 6) ~ 0.6 in both the cases. 



equal to those in the bulk. Eq. [7J is the desired van 
der Waals equation of state for a DLG, developed in this 
paper. When c+ = 0, i.e., ignoring nearest-neighbor cor- 
relations, Eq. [7] reduces to the usual MF expression of \i 
for an ELG []jf . 

The van der Waals equation of state notably does not 
have any free parameter. To demonstrate that Eq. [Jj 
indeed explains various features of DLGs, we assume a 
physically motivated approximate form of c+(n, E) where 
we use c+(n, 0) > c + (n 7 E), for any n ^ 0,1, as sub- 
stantiated in Fig. [3J This relation is expected on the 
ground that the driving field acts as an extra noise to 
break nearest- neighbor bonds [14[. For sufficiently small 
\K\, fj,(n,E) is a monotonically increasing function of n 
(the condition for non-monotonicity is discussed later). 
For n < 1/2, since c+(n, 0) > c + (n,E), the second term 
in the square bracket of Eq. [5] is positive and greater 
in equilibrium than in nonequilibrium. Since K < 0, 
we therefore get /i(n, E) > /j,(n, 0) for a given n. It 
implies that, if a driven system 1 with density m is 
in contact with the corresponding non-driven equilib- 
rium system 2 with density n-i where ni,ri2 < 1/2, the 
steady-state densities will be such that n\ < n-i- For 
n.1,712 > 1/2, this would be exactly the opposite, i.e., 
n\ > ri2- For n\= n<i = 1/2, equalization of the chemi- 
cal potential fi(l/2,E) = /i(l/2,0) implies that densities 
of a driven and the corresponding non-driven system in 
contact would indeed be the same. These results, which 
are derived above using the concept of equalization of a 
chemical potential for driven systems, are now verified in 
simulations presented in Fig. [3]where we plot the density 
n i of a driven system in contact with the corresponding 
non-driven system with density ri2 for various K. Note 
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FIG. 3: Density m of a 20 x 20 driven system as a function 
of density ni of the 250 x 250 equilibrium reservoir for vari- 
ous interaction strengths (particle exchange through a 2 x 2 
contact region). The simulations (points) are compared with 
the MF theory (lines) which is obtained using a particular 
choice of c+(n, 0) as given in Eq. 1101 and A(E) = 0.6, which 
corresponds to E — 6, as estimated in Fig. [2] Note that, in 
the simulations for K < K c ~ —1.30 as well as in the MF 
theory for K < J\^ IF ~ —1.22, there is a jump in the density 
rii when the reservoir density reaches ni = 1/2. 



that they are expectedly consistent with the particle-hole 
symmetry. 

Phase transition. - Now we describe the phase transi- 
tion, observed in the simulations, in terms of the thermo- 
dynamic potential /i, analogous to the description of the 
phase transitions in terms of the free energy in equilib- 
rium. As mentioned before, for sufficiently large \K\, fi 
can be a non-monotonic function of n. This is unphysical 
as one expects that, with increase in the density 712 (or 
equivalently /12) of the equilibrium system 2, the density 
rii should also increase. Here, the non-monotonicity of fi 
is the signature of the phase transition occurring below a 
critical value of K < K^ F . This is verified in the simula- 
tions in Fig. [3] where we plot the density n\ of a strongly 
driven system 1, with E\ = 6, as a function of density 712 
of the non-driven system 2. For K < K c ~ —1.3, there 
is a jump in the density m. The size of the jump goes 
to zero continuously as K tends to K c from below. The 
criticality condition in the MF approximation is given by 
(d/i/dn) = 0, i.e., 
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For any specific form of c+(n,E), the critical K^r can 
be found as a solution of K from Eq. [9] Generically, 
at K = K^ F and n = 1/2, (dM/dn) has a minimum 
touching the minimum in the function l/{n(l — n)}, i.e., 
the l.h.s. of Eq. [9] Therefore, writing [dM/dn} n — 1/2 = 
2[1— 4c + (l/2, E)], we obtain from Eq. [H]the critical inter- 
action strength K™ F (E) = -1/[1 - 4c+(l/2, E)]. Now, 




FIG. 4: The van der Waals equation of state: Chemical po- 
tential (j, as & function of density n. The curve for \i develops 
a kink for K < K™ F (E = 6) ~ -1.22. The size of the jump 
in density in the coexistence phase is obtained by using the 
Maxwell's construction (see the black curve). We have taken 
c(n, E) — A(E)c(n,0) (see Eq. [TOland the text below) with 
A(E) — 0.6 which corresponds to E — 6 (see Fig. [2]). 



using c+(n,0) > c+(n,E), we get K™ h (E) > Kf* (0) 
which has been observed in our simulations here as well as 
in simulations in the past [9|,|l0j|. Since the equalization of 
an intensive variable can be related to the existence of a 
generalized free energy f(n, E) [y, LL2j, we get n{n, E) = 
df/dn and therefore f(n,E) = /" n(n',E)dn'. When 
fj, is a non-monotonic function of n, f(n,E) would not 
be concave. But concavity of f(n,E) could be restored 
by the usual Maxwell's construction and the jump Sn 
in the density can be determined accordingly (see Fig. 
[4]). Since c+{n,E) is symmetric around n = 1/2 due to 
the particle-hole symmetry, one can see from Eqs. [7] and 
[8] that n(n, E) is anti-symmetric around n = 1/2. So 
the line // = /i(n = 1/2, E) gives equal areas which are 
bounded by the line and the fi(n, E) curve. The two den- 
sity values, where the line // = /i(n = 1/2, E) intersects 
the curve n(n, E) at the left and the right, correspond to 
the densities in the liquid and gas phase, respectively. 

We illustrate the above analysis by taking an approxi- 
mate form of the equilibrium correlation function 



c+(n, 0) ~ n(l — n) + 
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This form can be obtained using approximation of a dy- 
namical mean- field theory [15j . Then we assume, for sim- 
plicity, c+(n,E) = A(E)c+(n,0) where < A(E) < 1 
with A(0) = 1 in equilibrium (see Fig. [2]). In this special 
case, using Eq. [7J the chemical potential /i is plotted as 
a function of density n for various values of K in Fig. [4] 
with a specific choice of A(E) = 0.6 as estimated in Fig. 
[21 Note that the chosen value of A(E) corresponds to 
a strongly driven system with E = 6 since the driving 
field is much larger than the interaction strengths. The 
kink in \i appears at K^f F (E = 6) ~ —1.22 which indi- 



cates the onset of the phase transition. We also obtain 
the equilibrium MF value K^ F (E = 0) ~ -1.62. Evi- 
dently, both values are quite near to the corresponding 
known values [10J, K c ~ —1.3 in strongly driven case and 
K c ~ —1.76 in equilibrium. In Fig. [31 we plot the den- 
sity ni of the driven system 1 as a function of the density 
«2 of the non-driven system 2 where simulations and the 
MF theory agree quite well, except near the transition 
region. Moreover, as another validation of the MF the- 
ory developed here, in the bottom panel of Fig. Q] we 
plot the integrated compressibility I x vs. chemical po- 
tential /i obtained from the MF theory in the fluid phase 
with K = — 1 (i.e., K < K^f F ) where driven and equilib- 
rium cases are compared. The MF results, without any 
fitting parameter, are again in good agreement with the 
corresponding simulation results. 

Summary. - In conclusion, we have given an ap- 
proximate thermodynamic theory which captures vari- 
ous properties of paradigmatic driven lattice gases with 
attractive nearest-neighbor interactions remarkably well 
and, in particular, gives a consistent description of the 
phase transition from the homogeneous fluid phase to the 
coexisting phases of liquid and gas. Essentially, we have 
obtained a formalism to calculate a part of a putative 
nonequilibrium free energy which governs the phase tran- 
sition observed in various simulations. It still remains to 
be seen whether and how the long-range correlations ob- 
served in these systems [lOl Il6j , but so far ignored in our 
analysis, affect this thermodynamic theory of the phase 
transition. 

Importantly, our study opens up the possibility of ther- 
modynamic characterization of driven systems, in gen- 
eral, as the theory can in principle be extended to these 
systems which exchange a conserved quantity upon con- 
tact. However, the challenge in such an extension actu- 
ally lies in choosing a suitable contact dynamics so that 
equalization of a thermodynamic variable occurs. 

We thank R. K. P. Zia and R. Ramspcrger for discus- 
sions. 
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